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Abstract. We investigate thermoelectric properties of correlated quantum dots and molecules, 
fl} | described by a single level Anderson model coupled to conduction electron leads, by using 

Wilson's numerical renormalization group method. In the Kondo regime, the thermopower, 

S(T), exhibits two sign changes, at temperatures T = T\ and T = Ti > T\. We find that Ti is 

of order the level width F and Ti > T p ~ Tk, where T v is the position of the Kondo induced 

peak in the thermopower and Tk is the Kondo scale. No sign change is found outside the Kondo 

regime, or, for weak correlations, making a sign change in S(T) a particularly sensitive signature 

i-pj , of strong correlations and Kondo physics. For molecules, we investigate the effect of screening by 

<~] ■ conduction electrons on the thermoelectric transport. We find that a large screening interaction 

Q ' enhances the figure of merit in the Kondo and mixed valence regimes. 

(N 
> 

£NJ I 1. Introduction 

Materials with high thermoelectric efficiency are currently under intense theoretical and 
experimental investigation, largely due to the prospect of applications, e.g., for the conversion 
of waste heat into electricity, or for refrigeration and on-chip cooling in microelectronics [TJ [SJ . 
Apart from applications, thermoelectric materials also serve as an interesting testing ground for 
new theoretical approaches to thermoelectric transport in solids [31 HI [5]. Coulomb correlations 
and the Kondo effect give rise to large thermoelectric coefficients in strongly correlated bulk 
materials [Q El 13 [9]. They are also important in nanoscale semiconducting devices and 
in molecular transistors, posing interesting theoretical and experimental challenges for the 
understanding of electrical and thermal transport in such systems [TOj [Til HI E3 H5]IT5]. 

In this paper we address the thermoelectric properties of two related systems: a nanoscale 
size quantum dot exhibiting the Kondo effect, which we describe in terms of a single level 
Anderson impurity model with two conduction electron leads at fixed chemical potentials, and, 
a molecular transistor described by an extension of the above model to include a local screening 
of the Coulomb interaction on the molecule by the electrons in the leads. The quantum dots 
or molecules that we consider can be tuned from the Kondo to the mixed valence and empty 
orbital regimes by a gate voltage [17] IE]. Very recent experiments on nanoscale quantum dots 
|10| are beginning to probe the effect of Kondo correlations on the thermopower. 

We use Wilson's numerical renormalization group (NRG) method [1Q|, [20] to calculate the 
linear conductance, thermopower and thermal conductance as a function of temperature and 
gate voltage. The technique gives reliable results for transport properties in all regimes of 



in 



S 



Kondo 



Mixed valence 




10 10" 10 

T/T 



T/r 



Figure 1. Conductance, G, 

thermopower, S, and, thermal 
conductance, K e , versus T/T, for 
U = AT = 0.08 and gate voltages, 
v 9 = (e d +U/2)/T, in the Kondo (a- 
c) and mixed valence (d-f) regimes 
(Go = 2e 2 /h). The v g are indicated 
in top panels and increase by 0.25. 
Arrows show trends with increasing 
v g . T p : position of low temperature 
peak in S. Temperatures T\p at 
which S(T) changes sign are shown 
in (g) versus v g . (h): gate voltage 
dependence of S at representative 
temperatures T relative to Xi(0) 
and T 2 (0). 



interest [1]. The recently introduced full density matrix (FDM) approach [21] (see also [22]) is 
used. This allows calculations of dynamical properties at all excitation energies uj relative to 
the temperature T, thereby simplifying the calculation of transport properties which require 
knowledge of excitations, u>, above and below the temperature. 



2. Model and transport properties 

We describe both the nanoscale quantum dot and the molecular transistor by the single level 
Anderson impurity model with two conduction electron leads 
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Here, e a ka is the kinetic energy of electrons with wavenumber k and spin a in lead a = (L, R), £& 
is the local level energy, U is the local Coulomb repulsion and t a is the tunnel matrix element of 
the local level to conduction electron states in lead a = (L, R). -^screening = ^dc( n d — l)( n o — 1)> 
with rid the local level occupancy and no the local occupancy of the lead electrons, represents 
the screening interaction for the case of the molecule. We assume symmetric couplings of the dot 
to both leads and use the full width at half maximum, T = 0.02, as the energy unit throughout. 
Thermoelectric transport is calculated for a situation in which a small external bias voltage, 
SV = Vl — Vr, and a small temperature gradient ST is applied between left and right leads 
following |11[ 112] . To linear order, the following expressions for the electrical conductance, 
G(T), the thermal conductance, K C (T), and the thermoelectric power, S(T), are obtained 
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where I n , n = 0, 1, 2 are the transport integrals 



In(T) = lj ' du U n T(uj)(-^). 
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Figure 2. Temperature depen- 
dence of the thermopower, S(T), 
for U = 4T = 0.08 in the Kondo, 
(a), mixed valence, (b), and empty 
orbital, (c), regimes for increasing 
values of the screening interaction 
Udc- The case of large screening 
Udc = 0.16 > U is qualitatively sim- 
ilar to a large local electron-phonon 
coupling and results in destruction 
of Kondo correlations at finite gate 
voltages. Screening has a negligible 
effect in the empty orbital regime, 
but affects significantly S(T) in the 
other regimes. 



Here, e denotes the magnitude of the electronic charge and h denotes Planck's constant. The 
transmission function 7~(w) = ttTA(oj), where A(u) is the d-level spectral density. 

3. Results for quantum dots: Udc = 

The temperature dependence of the electrical conductance, G(T), thermopower, S(T), and 
electronic part of the thermal conductance, K e (T), are shown in Fig. [lk-c in the Kondo regime 
and in Fig. QJl-f in the mixed valence regime |23j. The thermopower exhibits two sign changes 



in the Kondo regime at T = Ti(v s ) and T = T2(v g ) and no sign changes on entering the 
mixed valence regime (see Fig. QJ). The gate voltage dependence of S in Fig. QJi shows three 
characteristic types of behaviour depending on the T relative to Ti(0) as explained in [23] . 

4. Results for molecules: Udc > 

For molecules in metallic break junctions, we expect additional effects to be important for 
thermoelectric properties, e.g. local phonon modes and screening by the conduction electrons 
in the leads. Here, we investigate the effect of the latter, and note that the effect of a strong 
local electron-phonon coupling can be approximately simulated by choosing a large screening 
interaction Udc > U. 

For moderate screening interactions Udc < U, we recover the expected downward 
renormalization of the local Coulomb repulsion U —> U* < U and an upward excitonic 
renormalization of the effective hybridization with increasing Udc leading to a weakening of 
Kondo correlations. Correspondingly, transport properties are significantly affected in the Kondo 
and mixed valence regimes, as shown in Fig. [2] for the thermopower and in Fig. [3] for the figure 
of merit. For large screening interactions Udc > U we find a destruction of Kondo correlations 
for gate voltages v g > away from particle-hole symmetry with an enhancement of the figure 
of merit for "Kondo" and "mixed valence" regimes (referred to the Udc = case) . 



5. Conclusions 

The thermoelectric power of, (i), quantum dots, and, (ii), molecules in break junctions, has been 
investigated using a single-level Anderson model. We showed that the thermopower exhibits sign 
changes as a function of temperature and gate voltage. These can be used as sensitive probes 
of Kondo correlations. Recent measurements on Kondo correlated quantum dots |10| are in 
qualitative agreement with our model calculations [23] . The effect of lead electron screening on 
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Figure 3. Temperature depen- 
dence of the "figure of merit" 
ZqT = GS 2 T/K C for U = AT = 
0.08 and increasing values of the 
screening interaction U^c i n differ- 
ent regimes (by reference to the 
Udc = case). The true figure of 
merit ZT = GS 2 T/(K C + K ph ) in- 
cludes a material specific phonon 
contribution, K p ^, which we have 
neglected, but needs to be included 
whenever phonons dominate ther- 
mal transport. Provided all scales 
are less than the Debye tempera- 
ture 6d, we expect that ZT « ZqT. 



molecular transport is shown to be significant in the Kondo and mixed valence regimes, where 
we find that an enhanced figure of merit can result. In future, it would be interesting to include 
a local electron-phonon coupling of the molecular state in order to investigate the effect of this 
interaction on molecular transport in detail [ 



Acknowledgments 

Support from the Forschungszentrum Jiilich (V.Z.), the SFB 767 (V.S.) and supercomputing 
support from the John von Neumann Institute for Computing (Jiilich) is acknowledged. 



References 

[I] Mahan G D 1998 Solid State Phys. 51 81 
[2] Kanatzidis M G 2010 Chem. Mater. 22 648 

[3] Zlatic V, Costi T A, Hewson A C and Coles B R 1994 Phys. Rev. B 48 16152 

[4] Costi T A, Hewson A C and Zlatic V 1994 J. Phys.: Condens. Matter 6 2519 

[5] Zlatic V and Monnier R 2005 Phys. Rev. B 71 165109 

[6] Arita R , Kuroki R K, Held K, Lukoyanov A V, Skornyakov S and Anisimov V I 2008 Phys. Rev. B 78 115121 

[7] Lackner R, Bauer E and Rogl P 2006 Physica B 378-390 835 

[8] Bentien A, Johnsen S, Madsen G K H, Iversen B B and Steglich F 2007 Europhys. Lett. 80 17008 

[9] Sales B C, Mandrus D and Williams, R K 1996 Science 272 1325 

[10] Scheibner R, Buhmann H, Reuter D, Kiselev M N and Molenkamp L W 2005 Phys. Rev. Lett. 95 176602 

[II] Kim T-S and Hershfield S 2002 Phys. Rev. Lett. 88 136601 
[12] Dong B and Lei X L 2002 J. Phys.: Condens. Matter 14 11747 

[13] Roch N, Florens S, Costi T A, Wernsdorfer W and Balestro F 2009 Phys. Rev. Lett. 103 197202 

[14] Parks J J, Champagne A R, Costi T A, Shum W W, Pasupathy A N, Neuscamman E, Flores- Torres S, 

Cornaglia P S, Aligia A A, Balseiro C A, Chan G K-L, Abruna H D and Ralph D C 2010 Science 328 1370 
[15] Koch J, von Oppen F, Oreg Y and Sela E 2004 Phys. Rev. B 70 195107 
[16] Reddy P, Jang S-Y, Segalman R A and Majumdar A 2007 Science 315 1568 
[17] Goldhaber-Gordon D, Gores J, Kastner M A, Shtrikman H, Mahalu D and Meirav U 1998 Phys. Rev. Lett. 

81 5225 
[18] Cronenwett S M, Osterkamp T H and Kouwenhoven L P 1998 Science 281 540 
[19] Wilson K G 1975 Rev. Mod. Phys. 47 773 

[20] Bulla R, Costi T A and Pruschke, T 2008 Rev. Mod. Phys. 80 395 
[21] Weichselbaum A and von Delft J 2007 Phys. Rev. Lett. 99 076402 
[22] Peters R, Pruschke T and Anders F B 2006 Phys. Rev. B 74 245114 
[23] Costi T A and Zlatic V 2010 Phys. Rev. B 81 235127 
[24] Leijnse M, Wegewijs M R and Flensberg K 2010 Phys. Rev. B 82 045412; Entin-Wohlman O, Imry Y and 

Aharony A 2010 Phys. Rev. B 82 115314 



